Contigs directed gene annotation (ConDiGA) for accurate protein sequence database construction in metaproteomics

Background Microbiota are closely associated with human health and disease. Metaproteomics can provide a direct means to identify microbial proteins in microbiota for compositional and functional characterization. However, in-depth and accurate metaproteomics is still limited due to the extreme complexity and high diversity of microbiota samples. It is generally recommended to use metagenomic data from the same samples to construct the protein sequence database for metaproteomic data analysis. Although different metagenomics-based database construction strategies have been developed, an optimization of gene taxonomic annotation has not been reported, which, however, is extremely important for accurate metaproteomic analysis. Results Herein, we proposed an accurate taxonomic annotation pipeline for genes from metagenomic data, namely contigs directed gene annotation (ConDiGA), and used the method to build a protein sequence database for metaproteomic analysis. We compared our pipeline (ConDiGA or MD3) with two other popular annotation pipelines (MD1 and MD2). In MD1, genes were directly annotated against the whole bacterial genome database; in MD2, contigs were annotated against the whole bacterial genome database and the taxonomic information of contigs was assigned to the genes; in MD3, the most confident species from the contigs annotation results were taken as reference to annotate genes. Annotation tools, including BLAST, Kaiju, and Kraken2, were compared. Based on a synthetic microbial community of 12 species, it was found that Kaiju with the MD3 pipeline outperformed the others in the construction of protein sequence database from metagenomic data. Similar performance was also observed with a fecal sample, as well as in silico mixed datasets of the simulated microbial community and the fecal sample. Conclusions Overall, we developed an optimized pipeline for gene taxonomic annotation to construct protein sequence databases. Our study can tackle the current taxonomic annotation reliability problem in metagenomics-derived protein sequence database and can promote the in-depth metaproteomic analysis of microbiome. The unique metagenomic and metaproteomic datasets of the 12 bacterial species are publicly available as a standard benchmarking sample for evaluating various analysis pipelines. The code of ConDiGA is open access at GitHub for the analysis of microbiota samples. Video Abstract Supplementary Information The online version contains supplementary material available at 10.1186/s40168-024-01775-3.

benchmarking sample for evaluating various analysis pipelines.The code of ConDiGA is open access at GitHub for the analysis of microbiota samples.

Background
The human body is composed of both human cells and many different microorganisms.Through host-microbiota interactions, microbes are closely associated with various diseases, including luminal diseases, immune diseases, metabolic diseases, and neurodegenerative diseases.To understand the role of microorganisms in human health and disease, it is necessary to characterize the changes in composition as well as the functional dynamics of microbiota.With the development of nextgeneration sequencing (NGS) techniques, metagenomics has greatly facilitated the study of microbiota [1,2].While metagenomics provides information on the composition and functional potential of microbiome, the method cannot reveal proteins actually expressed in the microbiome [3,4].Proteins, as biomolecules performing various functions within organisms, should be identified and quantified directly for microbiota function study.
During the past years, mass spectrometry (MS)based metaproteomics has been emerging as a powerful approach to understanding the functions of microbial communities [5][6][7][8][9].Compared to traditional proteomics of a single organism or simple mixtures, metaproteomics of microbiota faces challenges of high complexity and taxonomic diversity, wherein dozens or even hundreds of species can present in a sample with wide dynamic range variations, making the proteomic analysis extremely difficult [3,10].In shotgun proteomics, peptide identification relies on matching tandem mass spectra with protein sequences in a database.A crucial step in metaproteomics is to choose a suitable protein sequence database.An incomplete database can result in the missing of some key proteins, while a database with too large search space can result in limited detection sensitivity and high false discovery rates [11][12][13].
To date, the database-building strategies in metaproteomics include mainly (i) refining public proteome sequence databases, e.g., NCBInr [14] and UniProtKB/ Swiss-Prot [15], using mass spectrometric data [16][17][18][19][20]; (ii) filtering public proteome sequence databases using taxonomic information by 16S rRNA sequencing [13,21,22]; and (iii) constructing sample-specific proteome sequence databases by whole-genome sequencing [10,12,13].Among the different strategies, sample-specific database construction by metagenomics sequencing has been considered the best choice and employed in many metaproteomic studies [12,[23][24][25][26], wherein the databases contain only protein sequences specific to samples and hence can offer the best-fit search space to explore the protein expression of a particular microbiota and to improve the overall proteome coverage.However, accurate taxonomic annotation of metagenomics data, which acts as a key role in the construction of protein sequence databases, remains a challenging task [10,26].Currently, there are three major taxonomic annotation strategies, i.e., read-level annotation, gene-level annotation, and contig-level annotation.Different approaches for taxonomic annotation of genes result in significantly divergent results in the downstream metaproteomic studies [12,24].There is, to date, no consensus on a robust and reliable approach for gene taxonomic annotation.

Optimization of taxonomic annotations of genes based on simulated microbial communities
In this work, we developed an accurate gene taxonomic annotation strategy, namely contigs directed gene annotation (ConDiGA, MD3, Fig. 1).Contigs assembled from metagenomic sequencing results were firstly annotated.Afterwards, the most confident species were selected, and a gene-level annotation was performed based on the reference genome of the selected species.Two commonly used annotation strategies were compared with the Con-DiGA, including the MD1 pipeline of annotating genes directly against the whole bacterial genome database and the MD2 pipeline of taking the contig taxonomic annotation information as the corresponding gene taxonomic information (Fig. 1).In all three pipelines, different annotation tools, including BLAST [27], Kaiju [28], and Kraken2 [29], were compared.Referred to the work by Tanca et al. [24], we also built a meta6FT database by six-frame translation of the contigs using MMseqs2 [30], as well as a metaPA database where the amino acid sequences of the predicted genes were annotated using the UniProtKB/TrEMBL database using Protein BLAST.
We mixed 12 known bacterial species to derive a synthetic microbial community (Supplementary Table 1) to benchmark the different protein sequence database construction strategies.We also included the Uniprot reference proteome at the species level (UP-S), genus level (UP-G) and family level (UP-F) for comparison.The detailed information in terms of the number of genes annotated and species annotated of the protein sequence databases built by the three pipelines are shown in Supplementary Data 1, Supplementary Figs. 1, 2 and 3, and Supplementary Table 2.The Uniprot reference proteomes of the 12 species are shown in Supplementary Table 3. From the perspective of database annotation, MD1_BLAST acquired the highest annotation rate of over 90%.However, MD1_BLAST also had a high error rate in annotation (Supplementary Table 2), and 122 species were identified by the MD1_BLAST.The MD1 pipeline with all three annotation tools led to a significant number of false annotations at the species level.MD1_ Kaiju identified 157 species, and MD1_Kraken2 identified 52 species.MD2 was more accurate than MD1 in taxonomic annotation at the species level.MD2_BLAST identified 77 species, MD2_Kaiju identified 121 species, and MD2_Kraken2 identified 33 species.The MD3 is most accurate in taxonomic annotation at the species level.Since the genes were annotated based on a refined set of species, MD3 not only had a higher recall on gene annotations but also was robust against mis-assemblies.All the annotation tools with MD3 led to only 12 species being identified.However, it should be noted that MD3 with BLAST had one mis-annotation, wherein Klebsiella variicola was annotated instead of Lactobacillus acidophilus.Two close relatives of Klebsiella variicola, i.e., Klebsiella aerogenes and Klebsiella pneumoniae, were included in the synthetic community, which can explain the mis-annotation by MD3 with BLAST.
We then compared the different protein sequence databases by the different annotation strategies on the analysis of the metaproteomic data of the synthetic microbial community.It should be noted that the protein sequences in all the databases by MD1, MD2, and MD3 were the same, and only the annotation results were different, because the different pipelines shared the same gene prediction strategy.The proteins without annotation information at species level were labeled as unknown in the databases.Protein group identification results are shown in Table 1 and Supplementary Data 2. The identification results using the UP-S were considered as the reference.Generally, the numbers of protein groups identified for the 12 target species using the MD3-based databases were significantly higher than those using the MD1and MD2-based databases.As for the annotation tools, BLAST performed better overall, but Kaiju-MD3 managed to become the most efficient strategy, which led to the identification of 13,537 protein groups for the 12 target species.In contrast, the UP-S only led to 12,604 protein groups identified for the 12 target species.Figure 2A and B show the numbers of protein groups and peptides identified for the 12 target species using the databases of Kaiju with the MD3 pipeline (shown as MG), UP-S, UP-G, UP-F, and Meta6FT.The results by the UP-G, UP-F, and Meta6FT databases were not comparable to those by the MD3-based databases.For further analysis, we compared the protein groups identified for each species using MD3_Kaiju, MD3_BLAST, MD3_Kraken2, and UP-S (Fig. 2C).The numbers of protein groups assigned to each species using the MD3-based databases were basically consistent with those using the UP-S except for K. pneumonia and Lactobacillus acidophilus.The number of protein groups identified for K. pneumonia using MD3_Kaiju and MD3_Kraken2 was higher than that using UP-S.UP-S database was constructed using the reference stains, which can be different from the strains in our synthetic microbial community.By metagenomic sequencing, it is possible to identify some proteoforms with sequence variations not included in UP-S.Based on the benchmarking with the synthetic microbial community of 12 species, it was demonstrated that the MD3 pipeline is most suitable for annotating genes in the construction of metagenomics-derived database for metaproteomic analysis, and the MD3 strategy plus Kaiju served the best among all the tested methods.

Benchmarking the performance of ConDiGA on fecal sample
To test the performance of the ConDiGA pipeline in real microbial communities, a stool sample from a healthy volunteer was collected and processed for metagenomic and metaproteomic analyses.Protein sequence databases were constructed from the metagenomic sequencing results using the MD1, MD2, and MD3 pipelines with the annotation tools of BLAST, Kaiju, and Kraken2.The detailed information in terms of the number of genes annotated and species identified of the protein sequence databases are shown in Supplementary Data 3 and Supplementary Table 4. From the perspective of taxonomic annotation, MD3_Kaiju showed the highest annotation rate of over 60%, followed by MD2_Kaiju and MD3_BLAST with 51.41% and 44.67%, respectively.As depicted in Fig. 3A and Supplementary Data 4, the numbers of identified protein groups with taxonomic annotation employing the MD3 pipeline with all three annotation tools tremendously exceeded the ones with the MD1 and MD2 pipeline.The best performance was again obtained with the MD3 pipeline and the Kaiju-based annotation, where > 90% of the identified protein groups were successfully annotated at the species level.When considering the species identified by the different pipelines, the MD3 showed the best consistency with the different annotation tools, i.e., BLAST, Kaiju, and Kraken2.With the MD3 pipeline, MD3_Kaiju identified the largest number of species (143 species) and covered the majority of the species identified by MD3_BLAST (71.3%, 57/80) and MD3_Kraken2 (64.4%, 47/73) (Fig. 3B).
To further analyze the stool sample, we conducted a comparative investigation of functional and taxonomic features at metagenomic (MG) and metaproteomic (MP) levels.Metagenomic and metaproteomic results based on the MD3_Kaiju were annotated using KEGG Orthology (KO) through the GhostKOALA website.Figure 4A shows the comparison of the KO functional annotation between the MG and MP in terms of metabolism, genetic information processing, and environmental information processing.Notably, the KO annotation results by the two meta-omics methodologies were mostly consistent.Carbohydrate metabolism, (protein forming) amino acid metabolism, and energy metabolism were the top three abundant metabolism categories according to MP, while MG suggested carbohydrate metabolism, metabolism of other amino acids, and energy metabolism as the top three.This result highlighted the divergence between metagenomics and metaproteomics in revealing functional potential of microbiota.
The relative taxonomic abundance at the metaproteomic level was computed by summing the quantitative information of all the identified peptides of each species, and the relative taxonomic abundance at the metagenomic level was calculated by counting the number of genes assigned to each species (Supplementary Data 5).The cladogram in Fig. 4B depicts the relative abundance discrepancies between MP and MG by calculating the log 2 MG/MP abundance ratio.The results showed that the abundance differences existed at different taxa levels.MP and MG abundance were generally consistent at order level and levels beyond order, while their distinguishment appeared frequently below order level.As an example, in the Bacteroidaceae family, the family-level MG/MP ratio was close to one, while the MG abundance was higher in Phocaeicola genera and the MP abundance was higher in Bacteroides genera.These results indicated that within the same family, the relative abundance of different species by metagenomics and metaproteomics can display heterogeneity.The relative abundance by MG is closely related to the cell copy of a species, while the relative abundance by MP shows the total protein amount of a species.The relative abundance discrepancies between MP and MG have also been reported by Tanca et al. [31].

Benchmarking on in silico mixed metagenomic and metaproteomic data
To further demonstrate the annotation sensitivity and accuracy of the ConDiGA strategy, we computationally mixed the data of the synthetic microbiota community with the data of the stool sample.The metagenomic data were mixed at the reads level while the metaproteomic data were mixed at the MS raw data level.Metagenomics-derived protein sequence databases were constructed using the MD1, MD2, and MD3 pipelines with the annotation tools of BLAST, Kaiju, and Kraken2.The detailed information in terms of the number of genes and species annotated in the protein sequence databases is shown in Supplementary data 6 and Supplementary Table 5. Considering the gene annotation ratio, Kaiju-MD3 remained the best choice.All the annotation strategies based on the MD3 pipeline successfully recovered the 12 species from the mixed metagenomic datasets, while the rankings of the 12 species in the databases considering the number of annotated genes for each species were varied by the different methods (Supplementary Table 6).
The metagenomics-derived protein sequence databases were then used to analyze the in silico mixed MS data of the synthetic community and the fecal sample.As shown in Fig. 5, Table 2, and Supplementary data 7, in general, the MD3 strategy performed the best and significantly better than the MD1 and MD2.For BLAST, although the numbers of protein groups identified for the 12 species were comparable among the different pipelines of MD1, MD2, and MD3.The numbers of protein groups identified from the other species, i.e., from the fecal sample, were significantly lower with the MD1 and MD2 pipeline compared to the MD3 pipeline.In the mixed data, the 12 species can be viewed as the high abundant ones.Such results indicated that BLAST with MD1 and MD2 pipelines in protein sequence database construction were less sensitive to low abundant species compared to the MD3 pipeline.Among all the methods, the database by MD3 with Kaiju performed the best.We also merged the protein sequence databases of the synthetic microbiota community and the stool sample (merged database) for further comparison, where both databases were built by the MD3 pipeline with the Kaiju annotation using the metagenomic data separately.As shown in Fig. 5 and Table 2, the protein sequence database built from the in silico mixed metagenome data provided a performance similar to the merged database.We have also compared the species identified from the mixed metaproteome data or from the fecal metaproteome data.As shown in Fig. 5B, more than 78% species identified from the fecal sample were recovered from the mixed metaproteome data using the protein sequence database built by the MD3 pipeline with different annotation tools (BLAST, Kraken2, and Kaiju), which further illustrated the robustness and stability of the MD3 pipeline.Among the three tools, Kaiju with MD3 showed the highest recovery (81.1%, 116/143) of species from the fecal sample.All these results demonstrated the effectiveness of the MD3 strategy in protein sequence database construction from metagenome data and again illustrated that the MD3 with Kaiju can provide the optimal taxonomic annotation.

Discussion
In this study, we compared three taxonomic annotation pipelines to construct metagenomics-based protein sequence databases and assessed their performance in metaproteomics analysis using a lab-assembled microbial mixture, a real-life stool sample, and the in silico mixed data of the two samples.To date, one widely used approach in constructing metagenomics-based protein sequence databases is based on annotation of predicted individual genes [12,24,26] (MD1 in this study).We found that this pipeline could result in a large number of annotated species not actually present in the sample due to inaccurate annotation on short genes.As annotations on long contigs are expected to be more reliable than those on short genes, we further used the second pipeline (MD2) for taxonomic annotation, where the taxonomy information of contigs was passed on to its genes.The risk of this pipeline is that reads from different species can be wrongly assembled to one contig due to the partial sequence similarity among species.As shown in Supplementary Fig. 4, there was indeed a significant portion of contigs with more than one species information.Therefore, we proposed the third annotation pipeline (MD3), named ConDiGA, which first selected the most confident species from the annotation results of contigs, and then performed gene-level annotation against these most confident species.We compared the MD1, MD2, and MD3 annotation strategies with BLAST, Kaiju, and Kraren2 as the annotation tools through the synthetic microbiota community of 12 species.The results generally revealed that the MD3 pipeline outperformed the MD1 and MD2 strategies not only in terms of the numbers of protein groups identified, but also in terms of the sensitivity and accuracy of annotations.As for the three annotation tools, BLAST performed better overall but Kaiju-MD3 managed to become the most efficient strategy.This circumstance remained consistent for the real-life stool sample where the annotation rate of MD3 was significantly higher than that of MD1 and MD2, which proved that the MD3 pipeline remained its advantages when it came to the analysis of complex real-life samples.
To further benchmark the ability to deal with complex samples, we computationally mixed the metagenome and metaproteome data of the synthetic microbiota community and the stool sample and performed annotation by the aforementioned three pipelines.Again, the MD3 pipeline provided the most efficient and robust results.All 12 species were successfully recovered by the MD3 pipeline with different annotation tools.The MD3 combined with Kaiju achieved the largest number of protein groups annotated to the 12 species as well as to the other species from the fecal sample.

Conclusions
In summary, we proposed an accurate taxonomic annotation pipeline based on metagenomic data for deep metaproteomics analysis, namely contigs directed gene annotation (ConDiGA), and demonstrated the importance of correct gene annotation in constructing protein sequence databases for metaproteomics.In our Fig. 5 A Bar plots of the numbers of protein groups identified from the 12 species of the synthetic microbial community, from the other species, and without taxonomic annotation information.The protein sequence database was generated from the mixed metagenome data of the synthetic microbial community and the fecal sample using the MD1, MD2, and MD3 pipelines with different annotation tools.The MS data were mixed from the synthetic microbial community and the fecal sample.A merged protein sequence database of the 12 species and the fecal sample was used as the reference.B Venn diagrams of species identified from the in silico mixed metaproteome dataset and the fecal metaproteome dataset.The protein sequence databases were generated from the mixed metagenome data of the synthetic microbial community and the fecal sample or from the fecal sample only using the MD3 pipeline with different annotation tools strategy of ConDiGA, the assembled contigs were firstly annotated and then the most confident species from the annotation results were selected by considering genome coverage as well as taxonomic abundance.The predicted genes were then aligned against the reference genomes of the selected species.As for the performance evaluation of different annotation strategies, we compared current annotation strategies with ConDiGA on a synthetic 12-species microbiota community, a real-life stool sample, and the in silico mixed data of the two samples, using three state-of-the-art annotation tools, i.e., BLAST, Kaiju, and Kraken2.We found that ConDiGA surpassed all other annotation strategies in multiple areas such as annotation coverage, annotation accuracy, and annotation of low-abundance species.Our optimized taxonomic annotation pipeline can tackle the current problem of annotation reliability in metagenomics-derived protein sequence database and can promote the development of metaproteomics.The metagenomic and the metaproteomic data of the 12 species benchmark sample, as well as the ConDiGA pipeline, are publicly available, which can be used for the evaluation of various data analysis pipelines as well as the analysis of microbiota samples.

Bacteria culture and sample collection
The 12 microbial strains (Supplementary Table 1) used in this study were purchased from American Type Culture Collection (ATCC) or China Center of Industrial Culture Collection (CICC).After culture in tryptic soy broth (TSB), brain heart infusion (BHI), or De Man, Rogosa, and Sharpe agar (MRS), the concentrations of bacterial cells were measured using the plate counting method, and the optical density of bacterial solution determined at the wavelength of 595 nm (OD 595 ).The bacterial cultures were then washed twice with phosphate-buffered saline (PBS) at pH 7.4 (Solarbio, Beijing, China) and stored at − 80 °C.The 12 bacterial strains were mixed with different cell numbers (Supplementary Table 1) to construct a synthetic microbial community.The fecal sample was collected from a healthy volunteer and stored at − 80 °C for the subsequent experiment.

Metaproteomic sample preparation
The synthetic microbial community mixture was resuspended in a lysis buffer containing protease inhibitor cocktail (EDTA-free, 1 ×), 20 mM Tris-HCl pH 8.8, 8 M urea, and 1% sodium dodecyl sulfate (SDS) and then sonicated (50 W, 20 Hz, 10 min) using an ultrasonic crusher in an ice bath.Next, the obtained solution was centrifuged (12,000 g, 10 min, 4 °C) to remove cell debris.Proteins in supernatant were quantified using a Pierce BCA assay Kit (Thermo Fisher Scientific, Waltham, MA, USA) prior to lyophilization.For proteolysis, 300 μg of proteins were dissolved in 300 μL 8 M urea and 6.1 μL 0.5 M Tris-(2-carboxyethyl) phosphine (TCEP), followed Table 2 Identification results of protein groups from the in silico mixed metaproteome data of the synthetic microbial community and the fecal sample using different databases The in silico mixed metagenome data were used to generated the protein sequence databases.A merged protein sequence database of the 12 species and the fecal sample was used as the reference by incubation for 1 h at 37 °C with shaking at 600 rpm.Following this, 18 μL iodoacetamide (IAA, 0.5 M) solution was added into the dissolved protein solution and incubated at room temperature for 45 min (in darkness).The proteins were precipitated using 1.5-mL pre-cooled acetone (− 20 °C) for 4 h and washed twice with the precooled acetone.After drying at room temperature, the proteins were redissolved in triethylammonium bicarbonate (TEAB) solution (200 μL, 0.1 M).For trypsin digestion, 6 μg of trypsin (Beijing Hualishi Technology Ltd, Beijing, China) was added into the protein solution and incubated for 16 h at 37 °C with shaking at 600 rpm.Finally, the obtained peptides were transferred to Mon-oSpin C18 column (Tokyo Japan GL Sciences Inc) for desalting.After the desalting, peptides were quantified by Pierce quantitative colorimetric peptide assay (Thermo Fisher Scientific, Waltham, MA, USA).
For human stool sample processing, differential centrifugation was used to enrich microbial cells according to a previous report [7,32].Specifically, 0.5-g stool sample was mixed with 20 mL PBS and shaken for 30 min at 25 °C and 100 rpm.Then, the suspension was first centrifuged at 500 g, 4 °C for 5 min to remove the large particles.Next, the supernatant was collected and subjected to high-speed centrifugation (12,000 g, 10 min) to collect the precipitates.The precipitates from the 500 g centrifugation were subjected to the differential centrifugation strategy again.Then, the two precipitates from the 12,000 g centrifugation were combined as the final collected microbes.For protein extraction, the collected microbial precipitates were processed by liquid nitrogen grinding.The milled powder was dissolved in 0.5 mL lysis solution (100 mM dithiothreitol (DTT), 2% SDS, and 20 mM Tris-HCl pH 8.8).After heating at 95 °C for 30 min, the dissolved protein solution was centrifugated at 12,000 g, 4 °C for 10 min to collect the supernatants.To remove SDS, pre-cooled acetone solution (− 20 °C) was added and incubated under − 20 °C for 4 h.The precipitated proteins were collected by centrifugation and then washed twice using pre-cooled acetone.The obtained proteins were dried at room temperature and dissolved in 0.5-mL lysis buffer.The protein quantification method and proteolysis procedure were the same as those for the synthetic microbial community aforementioned.

LC-MS/MS analysis
For each sample, 10 μg of peptides was redissolved in 30 μL solvent A (0.1% formic acid in water) and analyzed by a nanoESI timsTOF pro mass spectrometer (Bruker, Bremen, Germany) with a nanoElute ® (Bruker Daltonics) LC system.For each injection of the timsTOF pro MS analysis, 200-ng peptides were separated by a C18 reversed phase analytical column with 1.6 μm resin (25 cm × 75 μm i.d., Ionopticks) by a 120-min gradient with phase A as 0.1% formic acid in water and phase B as 0.1% formic acid in 99.9% acetonitrile (Supplementary Table 7).The column flow rate was maintained at 300 nL/min with a column temperature of 50 °C.The instrument was operated in the data-dependent acquisition-parallel accumulation serial fragmentation (DDA-PASEF) mode with 10 PASEF scans per topN acquisition cycle and accumulation ramp times of 100 ms each.MS and MS/MS spectra were recorded from 100 to 1700 m/z with an ion mobility range (1/ K0) of 0.7-1.3versus/cm 2 .Charge was set to 0-5.The "target value" was set to 10,000.The dynamic exclusion was activated and set to 0.4 min.The quadrupole isolation width was set to 2 Th for m/z < 700 and 3 Th for m/z > 700.

Metagenomic sequencing and data analysis
DNA was extracted from the synthetic microbial community or the stool sample using HiPure Bacterial DNA Kits (Magen, Guangzhou, China) according to the manufacturer's instructions.The quality of the extracted DNA was detected using Qubit (Thermo Fisher Scientific, Waltham, MA) and Nanodrop (Thermo Fisher Scientific, Waltham, MA).The extracted genomic DNA was firstly fragmented by sonication to a size of around 350 bp, then end-repaired, A-tailed, and adaptor ligated using NEBNext ΜLtra DNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA) according to the manufacturer's preparation protocol.DNA fragments with a length of 300-400 bp were enriched by PCR.The PCR products were purified using an AMPure XP system (Beckman Coulter, Brea, CA, USA).The libraries were analyzed for size distribution by a 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA) and quantified using real-time PCR.Genome sequencing was performed on an Illumina Novaseq 6000 sequencer (Illumina, Inc., San Diego, CA, USA) using the pair-end technology.
Raw reads from the Illumina platform were filtered using FASTP (version 0.18.0) with the following criteria: (1) removing reads with ≥ 10% unidentified nucleotides; (2) removing reads with ≥ 50% bases having phred quality scores ≤ 20; (3) removing reads aligned to the barcode adapter.After filtering, resulted clean reads were assembled using MEGAHIT (version 1.2.9) [33] (with the parameters k-min = 21 and k-max = 141).Genes were predicted based on the final assembled contigs using Meta-GeneMark (version 3.38) [34] with the default parameters.The same method was applied to the synthetic microbial community dataset and the stool sample dataset.
In MD2, the assembled contigs were annotated instead of predicted genes using BLAST, Kraken2, or Kaiju with the same parameters as for MD1.Then, genes were annotated according to the species label of the contigs they belong to.
In MD3, the assembled contigs were annotated using BLAST, Kraken2, or Kaiju with the same parameters as for MD1.For the synthetic microbial community, the species with relative sequence abundance > 0.5% and coverage ≥ 0.1% were selected as the most confident species.For the fecal sample and for the in silico mixed data, the species with a relative sequence abundance ≥ 0.01% and genome coverage ≥ 0.1% were selected as the most confident species.Then, the predicted genes were annotated to one of the selected most confident species based on the best alignment from Minimap2 [35].
To build the Meta-PA database, we first predict the genes from the contigs using MetaGeneMark.The amino acid sequences of the predicted genes were annotated using the UniProtKB/TrEMBL database (downloaded on 15th of December 2022 from https:// www.unipr ot.org/ help/ downl oads) using Protein BLAST (BLASTP version 2.13.0 +).
To build the Meta-6FT dataset, we first processed the assembled contigs in an alternative way based on naïve six-frame translation.The Sequence Processor and Translator script from https:// cgpdb.ucdav is.edu/ DNA_ SixFr ames_ Trans lation/ was used to obtain the naïve six-frame translation.The command used included the following: python seqs_processor_and_translator_bin_ V118_AGCT.pyfinal.contigs.fafinal.contigs.xoutDNA 6 1 BIN 24.Results from all the six frames were combined to form one file meta.6FT.faa.Then, we used MMseqs2 [30] (version 13.45111) to align sequences to the Uni-ProtKB/TrEMBL database.The command used included the following: mmseqs createdb meta.6FT.faatarget_seq; mmseqs search target_seq TrEMBL_mmseqs/uniprot_ trembl_db results./tmp/;mmseqs createtsv target_seq TrEMBL_mmseqs/uniprot_trembl_db results meta.6FT_annot.tsv-full-header.Each gene was annotated with the species having the longest alignment with over 10% alignment fraction and over 10% sequence identity.

Metaproteome data analysis and bioinformatic analysis
Uniprot-based databases were built from the reference proteome of the used microbes at the species level, genus level, or family level.For species level Uniprotbased database, the reference stains for each species shown in Supplementary Table 3 were used to construct the database, while for the genus-level or familylevel Uniprot-based databases, all strains included in the corresponding genus or family were used to construct the databases.The metagenomics-based databases were translated from the annotated genes using MetaGeneMark (version 3.38) [34].Database searching for all metaproteomic data was carried out using PEAKS Studio (version X pro, Bioinformatics Solutions Inc., Canada).The metaproteomic data were searched with the following parameters: precursor ion tolerance 15 ppm, fragment ion tolerance 0.05 Da, maximum of 2 missed cleavage sites, carbamidomethylation (C) of cysteine as fixed modification, oxidation of methionine and deamidation (NQ) as variable modification, trypsin as proteolytic enzyme, and 1% false discovery rate threshold at both peptide and protein group level.KEGG annotation was performed using the Ghost-KOALA web server.Data visualization was conducted with R (version 3.5.1,https:// www.r-proje ct.org/) using the packages of ggplot2 (https:// github.com/ tidyv erse/ ggplo t2) and AntV G2 (https:// github.com/ antvis/ g2).

Fig. 1 A
Fig. 1 A Benchmarking different protein sequence databases from metagenomics or from Uniprot reference proteome using a synthetic microbial community of 12 species.B Pipeline of three annotation strategies (MD1, MD2, and MD3)

Fig. 2
Fig.2Numbers of A protein groups and B peptides identified for the 12 target species from the synthetic microbial community using different databases.C Numbers of protein groups identified from each species using the databases constructed with the MD3 pipeline and the UP-S reference database

Fig. 3 A
Fig. 3 A Bar plots of the numbers of protein groups with and without taxonomic annotation identified from the stool sample using the MD1, MD2, and MD3 pipeline-based databases.B Venn diagrams of the species identified from the stool sample using the MD1, MD2, and MD3 pipelines-based databases

Table 1
Identification results of protein groups for the 12 target species from the synthetic microbial community using different databases Total 13,852 13,852 13,852 13,852 13,852 13,852 13,852 13,852 13,852 13,852 3784 12,604 5999 6004